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Abstract 

Scalar fields appear in many theories beyond the Standard Model of particle physics. 
In the early universe, they are exposed to extreme conditions, including high temper¬ 
ature and rapid cosmic expansion. Understanding their behavior in this environment 
is crucial to understand the implications for cosmology. We calculate the hnite tem¬ 
perature effective action for the field expectation value in two particularly important 
cases, for damped oscillations near the ground state and for scalar fields with a flat 
potential. We find that the behavior in both cases can in good approximation be 
described by a complex valued effective potential that yields Markovian equations 
of motion. Near the potential minimum, we recover the solution to the well-known 
Langevin equation. For large held values we hnd a very different behavior, and our 
result for the damping coefficient differs from the expressions frequently used in the 
literature. We illustrate our results in a simple scalar model, for which we give an¬ 
alytic approximations for the effective potential and damping coefficient. We also 
provide various expressions for loop integrals at hnite temperature that are useful for 
future calculations in other models. 
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1 Introduction 

1.1 Scalar fields and the early universe 

Scalar fields play important roles in particle physics and cosmology. They appear in many 
theories beyond the Standard Model of particle physics (BSM). For instance, the low- 
energy effective description of string theory and other theories involving extra dimensions 
generally include numerous scalar helds, called moduli, that parameterize the properties 
of the compactihed internal dimensions (see e.g. DEI for review). Scalar helds can drive 
cosmic inhation (inhaton) |3H6], explain the strong CP-problem (axion) [7j, induce phase 
transitions [HHn], be Dark Matter candidates [T2tiT5] or Dark Energy [TGHIS]. The recent 
discovery of a Higgs boson [201 [21] conhrms that they indeed exist in nature. 

In cosmology involving scalar helds, it is crucial to understand their time evolution in 
the early universe, where they are exposed to extreme conditions including high tempera¬ 
ture, large energy density and rapid cosmic expansion. Their evolution in a time-dependent 
background provided by the primordial plasma and the cosmic expansion is a nonequilib¬ 
rium process. The initial conditions for this process depend on the model with the specihc 
scalar helds, and the mechanism that set the initial conditions is unknown in almost all 
models. Many observational hints point towards cosmic inhation, but alternative models 
such as bouncing universe scenarios (see [22l [23] for review) have been also studied, though 
there are still controversies [SUES]. The bouncing models also invoke scalar helds, typi¬ 
cally inspired by string theory, to resolve the cosmic singularity by a bouncing phase in 
the very early universe [26l [27] . In general, the held value at initial time is far away from 
the potential minimum. If the ehective potential is sufficiently steep, the held will relax to 
its minimum very quickly. However, there are numerous examples for scalar helds with a 
hat potential that evolve slowly compared to the time scale related to the propagation and 
interactions of individual particles, including the moduli, inhaton and axions. With the 
present work, focusing on the hnite-temperature ehects in a thermal bath, we aim to make 
progress towards a quantitative understanding of the nonequilibrium dynamics of scalar 
helds in the nontrivial backgrounds of the early universe. 

1.2 Motivations and assumptions 

Let us consider a real scalar quantum held 0 with a symmetry under 0 —)■ —0 that couples 
to other degrees of freedom, which we collectively refer to as At this stage we do 
not make any assumptions about the spin or interactions of the helds A*. If we restrict 
ourselves to renormalizable operators, then the most general Lagrangian reads 



0^00^0 - H(0) - + £, 


( 1 ) 


with 


V(4>) 


m 
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Here 0[Xi] is a sum of operators that depend on combinations of the helds dfj only and 
is a Lagrangian that specihes their masses and interactions amongst each other. The 
generalization to non-renormalizable operators is straightforward. If we neglect the inter¬ 
actions with Xi for a moment, then the zero mode of 0 in a Friedmann-Robertson-Walker 
spacetime would classically follow the equation of motion 

^ + 3H4> + = 0, (3) 

where H is the Hubble constant. The full quantum held 0 can be decomposed into 

(j) = (p + r]. (4) 

Here (p is the expectation value of the held 0, while the held rj contains quantum huc- 
tuations, such as single particle excitations. If we switch on the interactions with the 
background medium and take into account radiative corrections, one could expect p to 
follow an equation of motion of the form 

(p {3H -\- T^p)ip dipV{p) = 0. (5) 

Here V{p) is an ehective potential for p that includes radiative and thermal corrections, 
while T^ is the rate at which energy is dissipated from p to other degrees of freedom. 
The equation (|^ looks physically intuitive, but should be rigorously justihed when being 
applied. It cannot always be valid because the nonequilibrium dynamics of interacting 
quantum helds involves non-Markovian memory kernels. It is the goal of this work to 
derive an ehective equation of motion of the form ([^ from hrst principles, and to determine 
V{p) and explicitly in an illustrative simple model. Together with the ehective kinetic 
equation for the //-propagators derived in [28], this allows to understand the dynamics of 
the scalar held in terms of Markovian equations. Our derivation of (|^ uses a somewhat 
diherent method than previous approaches in [29II85] . but yields a result that is consistent 
with previous works in the limit of small held values, where a comparison can be made. 
The explicit expression for at large held values, however, is considerably diherent from 
those found in the literature. 

In order to obtain a Markovian equation of motion for (p, we have to make several 
approximations. As usual in the derivation of ehective kinetic equations [36] , these rely on 
a separation of time scales. In closed systems, such a separation is automatically given if 
one assumes a sufficiently weak coupling. In this case the rates at which particles scatter 
are much smaller than their ehective masses. This implies that macroscopic properties of 
the system change much slower than the time scale associated with individual scatterings. 
It also ensures that the system can be described as a dilute gas of well-dehned quasiparticles 
with ehective masses M^, Mx^ and widths T^, VXi with T,, M^, Vx^ ^ Here is 
the pole mass of the resummed //-propagator (the same for and in general ehective 
mass and width both depend on momentumj^ 

^In thermal equilibrium the dispersion relations and widths of quasiparticles can be read off from the 
pole structure of the real time propagators m- This definition can be generalized to situations far from 
equilibrium, see e.g. [28]. The propagators in a medium can in principle have a complicated pole structure 
due to the appearance of collective excitations, such as luons [38]. We ignore this subtlety here. 
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In the early universe, Hubble expansion brings in an additional time scale. In the 
present work we make two basic assumptions 

1) The effective masses of all particles that couple to (p are larger than the rate of Hubble 
expansion, i.e. H < 

2) The interactions in the primordial plasma are sufficiently weak to guarantee T,,, T^ -C 
M^, (and analogous for all other helds ^i). Here is the local curvature of V((p). 

3) The change in the effective masses and widths on microscopic times scales ~ l/Tf,?, 

should be small compared to these quantities themselves, i.e. M^, T^ T^M,, 

etc. 

For practical purposes we make two additional assumptions that are not required for the 
derivation of a Markovian equation of motion, but lead to considerable simplihcations 

4) 0 is the only held with a non-zero expectation value, i.e. (Ai) = 0. 

5) The thermodynamical state of the constituents Aj of the primordial plasma can be 
characterized by a single parameter, an effective temperature T. This is usually 
justihed if r;^’. > H. 

Assumption 1) allows to use Minkowski-space propagators in loop diagrams. Although H 
is absolutely crucial in the kinetic equation for (p on macroscopic time scales, it can be 
neglected when determining the transport coefficients in this equation from microphysics. 
Without this assumption the computation of loops is technically very difficult [2SH13]- 
Physically this assumption can be interpreted in terms of a separation of time scales: The 
microphysical time scales in processes described by loop diagrams (which are related to the 
inverse of the masses) must be much faster than the macroscopic Hubble expansion, such 
that the background metric is approximately constant during the duration of individual 
processes. Assumption 2) is simply necessary for any perturbative treatment to be valid. It 
also ensures that the effective masses (which include contributions that depend on ip and T, 
which evolve with time) are approximately constant on the microscopic time scale that is 
relevant in loop integrals. Assumption 3) is the physical “adiabaticity assmption” that the 
properties of (quasi)particles do not change signihcantly within the duration of individual 
microscopic processes. Assumption is for simplicity only, and it is straightforward to 
relax it. Assumption 5) implies that Aj have reached some degree of kinetic equilibration. 
Though thermalization in general is a complicated problem [441 - 154] , this assumption seems 
at least qualitatively justihed in the case under investigation here. If 0 has a hat ehective 
potential, then this means that it must have rather weak interactions (weaker than the 
interactions of the Aj amongst each other). This implies that the time scales l/T^, 1/if on 

^This assumption is usually not fulfilled during inflation in the most popular models. It can, however, 
apply in several other applications that we have in mind, such as the late time behaviour of moduli, warm 
inflation scenarios and, more general, the behaviour of an order parameter in a hot plasma that is not in 
de Sitter space. 
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which if evolves are much longer than the time scales l/F Xi on which the Xi relax to local 
thermal equilibrium. Making the same assumption for rj is far less justihed. For simplicity, 
we will nevertheless assume that the occupation numbers of rj modes can be characterized 
by the same effective temperature as the rest of the plasma. This assumption, which has 
been made in almost all past studies (often without mentioning), allows to use equilibrium 
propagators for all helds. 

The setup dehned by the assumptions l)-5) is not sufficient for a complete understand¬ 
ing of scalar helds in the early universe. In particular assumptions 1) and 5) considerably 
constrain its applicability. Moreover, the properties of the primordial plasma in realistic 
models of particle physics is far more complicated than the toy model presented in section]^ 
However, the analytic expressions we hnd are derived more systematically and consistently 
than any comparable results in the literature we are aware of. Moreover, in the course 
of their derivation we hnd analytic estimates for various integrals in hnite-temperature 
held theory that will be very useful for calculations in more realistic models. Finally, in 
spite of the comments above, there are at least two cosmological problems in which the 
assumptions l)-5) are widely believed to be justihed: warm inhation [33] and the fate of 
moduli at late times [33] . 

1.3 Overview of this article 

This paper is organized as follows. In section we briehy review the elements of nonequi¬ 
librium quantum held theory to set up the theoretical framework and notations for the 
following sections. In section we derive the approximate ehective action and the equa¬ 
tion of motion for the expectation value <^9 in a thermal bath. These results are illustrated 
in a simple scalar model in section and we present approximate analytic estimates for 
ehective potentials and dissipation coefficients. In section we summarize the results of 
this work and discuss some related issues in comparison with previous studies. Section 
is devoted to conclusions and outlook. In the appendices we calculate the spectral self¬ 
energies from setting-sun diagrams, which are needed for the computation of the dissipation 
coefficients in section HI 

2 Elements of nonequilibrium quantum field theory 

The standard methods to calculate S-matrix elements in particle physics are not suitable 
to describe systems far from equilibrium at large density because there is no well-dehned 
notion of asymptotic states, the properties of elementary excitations may signihcantly diher 
from those of particles in vacuum and classical particle number in general is not a suitable 
quantity to characterize the system. 

However, all observables can be expressed in terms of time-dependent correlation func¬ 
tions of the quantum helds without reference to asymptotic states or free particles. In 
most practical applications, all relevant information is contained in the one- and two-point 
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functions. The expectation valne (p can be identified with the one-point fnnction, 


ip = {4>). 


( 6 ) 


That is, the average 


(...) = Tr(p...) 


(7) 


is taken over qnantnm as well as statistical flnctnations encoded in the density operator 
Q. We will in the following always assnme that (p = 0 in the thermodynamic gronnd state 
and stndy how the system relaxes to this state if (p 7 ^ 0 at initial time. There are two 
independent two-point fnnctions. Common choices for these are the connected Wightman- 
fnnctions 


A>(a;i,a;2) = {(t){xi)(t){x2)) - ^{xi)ip{x2), (8) 

/\^{xi,X2) = {(t){x2)(t){xi)) - (p{xx)(p{x2) ( 9 ) 

or their linear combinations 

^~{xi,X2) = i {A>{xi,X2) - A<{xi,X2)) , A+(a;i,a;2) = ^ (A>(a;i, 0 : 2 ) + A<(a;i,X 2 )) .(10) 

These have a clear physical interpretation. The Fonrier transform 

pr,{p-,Xi +X 2 ) = j ^^(^1 “^ 2 ) A~{xi,X2) ( 11 ) 

of A“ is the spectral density pr^, the poles of which in po determine the spectrnm of qnasi- 
particles in a medinm. The Fonrier transform of A’*', on the other hand, characterizes the 
occnpation nnmbers of different held modes and can be interpreted as a held theoretical 
generalization of the classical particle distribntion fnnction, see e.g. discnssion in |2H]. In 
most practical applications, <p, A+ and A“ contain all relevant informations. 

2.1 The Closed-Time-Path formalism 

A convenient framework to obtain eqnations of motion for p, A+ and A“ is ohered by 
the Closed-Time-Path (CTP) formalism of noneqnilibrinm qnantnm held theory [53-1601 . 
In this formalism, correlation fnnctions are dehned with time argnments on the contonr 
shown in hgnre which starts from an initial time ti + ie, rnns parallel to the real axis 
to tf + ie, tnrns aronnd to tf — ie and retnrns to t^ — ie. We consider the limit —)• 00 , 
e —)• 0 and ti —)■ — 00 , where bonndary conditions at hnite time to can be imposed by 
external sonrces localized at to |2H]- In the following we dehne the qnantities reqnired for 
the present analysis. A more detailed introdnction to the CTP-formalism can be fonnd in 

[Ml Eg. 

The standard strategy to calcnlate correlation fnnctions is to split the held ef) on the 
contonr into (j)'^ and (j)~, which denote the held with time argnment on the “forward” and 
“backward” part of C, see e.g. [62]. One can dehne a Feynman propagator (Tc 0 (a:i) 0 (a: 2 )), 
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Figure 1: The contour C in the complex time plane. 


where Tc indicates time ordering along the contour path C. It can be decomposed into the 
correlation functions for these fields as 


Ai^^{xi,X2) = {T(j)'^{xi)4>^{x2)) - (p{xi)ip{x2), ( 12 ) 

^^~iXuX2) = {(p~{x2)(p'^ixi)) - ip{xi)ip{x2), (13) 

A-^{xi,X 2 ) = {(j)~{xi)(j)'^{x 2 )) - (p{xi)(p{x 2 ), (14) 

A"“(a;i,a;2) = {'f(j)~{xi)(l)~{x2))-ip{xi)ip{x2), (15) 


where T is the usual time-ordering and T is anti-time-ordering (because the backwards part 
of the contour runs from positive to negative infinity). Considering that time arguments 
on the forward branch are always “earlier” along the contour than those on the backward 
branch, one can easily identify 


A+ (a;i,a;2) = A<(xi,a:2) , A +(xi, 0 : 2 ) = A>(a;i, 0 : 2 ). (16) 


The same decomposition can be performed for the self-energies, see appendix of [63]. Then 
the combination (pc = (0^ + 0~)/2 is identified with a physical field and 0 a = 0 "'' — P~ is 
treated as a response field. After integrating out the bath fields A) and 0 a, one obtains an 
effective action for 0^, from which one can obtain its equation of motion. It can be expanded 
in powers of pc- At lowest non-trivial order one finds the Langevin type equation (46), 
including higher powers of pc in the effective action leads to nonlinear interactions between 
the different held modes (“multiplicative noise”), see e.g. jail Eg. By taking quantum and 
thermal averages over products of 0c, one can then calculate all correlation functions for 


2.2 Equation of motion for the background field ip 

In the present context we are, however, mainly interested in the one-point function p. 
We can therefore simplify the procedure and directly write down an effective action for ip 
(instead of p) on the closed time path C (MHSTj, 


n=l 

Here the label c indicates that the time integral is to be taken along C. The (xi,..., Xn) 
are the amputated vertex functionals. The expansion in ([T^ is around the minimum at 



/ d'^xi... —..., x„)(p(xi)... (^(x„). 

Ic nl 


(17) 






9 ? = 0 , and the coefficients are taken at (^ = 0 in the sense that they contain the held 
7] dehned in Q. However, for a consistent expansion in the small conplings reqnires the 
use of full (resumed) propagators in the loops because the T- and 93 -dependent corrections 
to effective masses bring powers of the couplings into the denominator of the propagators 
(similar to the situation of gauge theories at high T [ 68 ]). This is done in the 2PI formalism 
of non-equilibrium quantum held theory, see e.g. [HI]. Hence, the implicitly depend 
on 93 via the ehective masses. 

With ([^ in mind, we expand (17) up to fourth power. F^^^ is eliminated by the sta- 
tionarity condition for vanishing external sources. 




(18) 


F^^) is simply the inverse dressed propagator = 5c{xi — X 2 ) + 

Ilc{.Xi,X 2 )- F*^^) vanishes due to the symmetry 0 -H- —0 in (|^ and is the amputated 
four-point function. Due to energy-momentum conservation we can replace 


d'^Xi... (7^X4 — F^^^(a;i, ..., X4)93(a;i)... 93 ( 0 : 4 ) 


-)■ 


d‘^Xi(fx2—^^^{xif[X^ + Iic{xi,X2)](p{x2y 


(19) 


i.e. express the ehective action in terms of an ehective potential V[ 93 ] = —( 5 *^^)( 0 )F[ 93 ]. The 
function H is, for example, given by the diagram in hgure 

We now perform a spatial Fourier transform to the mixed representation with arguments 
it,p) and decompose 93 into modes p, which in general are coupled to each other by the 
integrals in ( [l7| . We are interested in the case, in which only p = 0 mode of 93 signihcantly 
dihers from zero, and we can thus perform the spatial integrals analytically. Using the 


stationarity condition (18) we can hnd an equation of motion for 93 on C, 


(dt + dt'Ilcit, t')(p{t') + dt'Ilcit, = 0 . ( 20 ) 

Here we have suppressed the momentum index and all functions are to be evaluated at 
p = 0. Note that this equation does not contain a Brownian noise term ^ on the RHS, as 
one might have expected. The reason is that this term comes from thermal fluctuations, 
and the expectation ([^ that is taken in the dehnition 0 of 93 includes an average over 
thermal fluctuations. A noise therm S, does appear in the (Langevin type) equation of 
motion (46) for 0 if one only traces over the bath degrees of freedom [291 - I54] . While these 


fluctuations cancel out in the equation for the one-point function 93 , they crucially affect 
n-point correlation functions of 0 with n > 1 . 


2.3 


Equations of motion for the fluctuations 77 and perturbation 
theory 


In order to calculate the coefficients H and H in (20) one has to evaluate loop integrals 
with Tj and Aj running in the loop. This requires knowledge of the propagators. Their 
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equations of motion are the Kadanoff-Baym equations (see (SH]), which can be derived in 
a similar manner as (20), see e.g. [33], 


r-ti 


(□i + (x,,i2) = - / / dt'Il^{xj,x')A^(x',X2) , 


( 21 ) 


'^2 


nil 


{ni + ml)A+{xi,X2) = - / dV / dt'U {xi,x')A+{x',X2) 




+ / dV / dt'Il+{xi,x')A {x',X2) , 


( 22 ) 


where mrj is a vacuum mass of //-quanta]^ and the self-energies are defined, equivalently to 
(10), as 

Il~{xi,X2) = Il>{xi,X2) - Il<{xi,X2) , n+(a:i,a:2) = -- (n>(a;i,a;2) -h Il<{xi,X2)) • (23) 
The usual retarded self-energy can be identihed with 


n^{xi,X2) = 9(ti - «2)n (a:i,X2). 


(24) 


Approximate Markovian equations - Correlation functions and self-energies for the 
bath helds A* can be dehned analogously. In general, the coupled second order integro- 
differential equations (21) and (22) are difficult to solve. However, in a weakly coupled 
theory, and in a background that changes slowly compared to the time scale of individual 
particle scatterings, the system can be described as a gas of quasiparticles. In a homoge¬ 
neous medium, the spectral density defined in ( pT| only depends on time, p,,(p; x) = Pr^{p; t). 
Let Hp(t) be a pole of Pr^{p;t), with 


= Refl,,(t) and = 2Imf2,j(f). 


(26) 


Both are in general time-dependent due to the interaction with the medium. In weakly 
coupled theories one observes the hierarchy 


r,(t) < n.it). (26) 

In Minkowski space this is enough to justify a “quasiparticle approximation”, i.e. to under¬ 
stand most properties of the system by picturing it as a gas of weakly coupled quasiparticles. 
In a time-dependent background we must make the additional assumption that the above 
quantities change adiabatically, leading to a more precise formulation of our assumption 

3), 

(irj and T^ (27) 

^Note that niri is in general different from m,p since it obtains y-dependence due to self-interaction of 
(/): For example, for self-interaction. 
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Under these assumptions Pr^{p',t) at fixed t features peaks of width ~ r^(t) at energies 
Po ~ which can be interpreted as quasiparticle-Yesonaiices,a.iad loop integrals are 

typically dominated by the regions near these poles. Here (and in the following) we omit 
the spatial momentum index. In this situation, one approximately finds |28] 

sin 

COS 

Here f{t) is a generalized distribution functions, which follows the Markovian equation of 
motion 





dj{t) = [i + f{t)]r<{t)-f{t)r>{t) 


with the “gain” and “loss” term^ 


(!>.*) 


2nn{t) 


P0=O,,(t) 


see also CD]. 


(30) 


(31) 


Special case - thermal equilibrium - If the background medium is in thermal equilib¬ 
rium, the functions A~ and H” are independent of the coordinate Xi + X 2 [33]. Moreover, 
the Fourier-transforms of the self-energies in the coordinate Xi — X 2 are related by the 
Kubo-Martin-Schwinger (KMS) relation 


n<(p) = e-wi’'n>(p), 


(32) 


nt(p) = -M 3 + feiro) n (p), 


(33) 


or equivalently 


where fsipo) = — 1) ^ is the Bose-Einstein distribution. In this case, pr, reads [33 

-2Imn^(p) -h 2poe 


PriiP) = 


{pI - m2 - p2 - Ren^(p))2 (Imn^(p) poe^ ’ 


(34) 


where e is an infinitesimal parameter. Note that the self-energies and all depend 
on T. It is clear from (34) that, if condition (26) is fulfilled, the quasiparticle dispersion 
relation (or “mass shell”) Qri is essentially fixed by Ren^(p) via the condition 


pj - - mi - Ren®(p) = 0, 


( 35 ) 


^Here we ignore subtleties about finite initial time in definitions. 
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while Imn^(p) gives the thermal width, 


r.= 




U ^'fj U 

which is just the difference of gain and loss term. Ren^(p) contains a zero temperature 
divergence that can be absorbed into the physical mass in the same way as in vacuum 
[321133] , In the following we interpret as the physical mass in vacuum, with Ren^(p) 
being hnite. The effective mass can be dehned as flri for vanishing spacial momentum 
p = 0. It is very useful to £t a Breit-Wigner function to (34) near the pole, which is 
parametrized by T^ and as 

Po^v 




(p) = 2Z 


{pl - Qiy + {poT, 


+pr(p)- 


(37) 


Here is the continuous part of Pr^{p). To obtain the correct residue, we introduced 

the parameter 


Z = 


1 - 


1 aRen^(p) 


2a 


dpo 


-1 


(38) 


J PO — ^r] 

For weak coupling, it is often sufficient to use the “zero width limit”, 


pj(p) = 27rZsign(po)^(Po “ + P““*(p)- (39) 

Other useful relations that follow directly from the dehnitions of the correlators are 

n"(p) = 2ilmU^{p), (40) 

^^{p) = fB{Po)pv{p) , A^(p) = [1 + fB{po)]pv{p), (41) 

and equivalently for the self-energies of other fields. Note that all the above functions 
depend on T, which we have not made explicit here to simplify the notation. Correlators 
and self-energies for all other helds are defined analogously. 


Perturbation theory - The perturbative expansion for correlation functions in a medium 
can be performed in terms of Feynman diagrams, analogous to the vacuum case. One differ¬ 
ence is that one has to use the propagators or, often more conveniently, their transforms 
in Wigner space cni. Knowledge of A^ allows to determine all other correlation functions, 
in particular 


A++(a;i,a;2) 

2 

= A+(a;i,a; 2 ) - -sign(ti - t 2 )A“(a;i, a; 2 ). 

(42) 

A--{xi,X2) 

= A+(a;i,a; 2 ) + -sign(ti - t 2 )A“(a:i, 0 : 2 ), 

(43) 

A“+(a;i,a;2) 

= A+(a;i,X2) --A-(a;i,X2), 

(44) 

A+“(a;i,a;2) 

= A+(a;i,a;2) +-A-(a:i,a;2). 

(46) 
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Another difference from the vacuum case is the fact that (j)'^ and (j)~ formally have to be 
treated as two different helds. There are no vertices that directly connect a 0"'“ line to a (f)~ 
line, i.e. the vertices are always either of “+” type or type, but the propagators (44) 
and (45) mix the fields and can connect a vertex to a “+” vertex and vice versa, see 


e.g, BII for a detailed discussion. 


3 Effective action and equation of motion for ip 


The known equations (28), (29) and 


allow to describe the dynamics of the held 
huctuations in terms of the Markovian equations, which can be physically interpreted as 
effective Boltzmann equations for quasiparticles in the plasma. In the following sections 
we derive an approximate Markovian equation for the held (p, so that the system can be 
completely described in terms of Markovian dynamics. 


3.1 Small field values: Brownian motion and Langevin dynamics 

We consider the case in which all degrees of freedom except (p are in thermal equilibrium. 
If (j) at initial time is very close to its thermodynamic ground state, then the ehect that 
this deviation has on the thermal bath (“back reaction”) is negligible. This in particular 
requires that ip be small enough that we can use a quadratic approximation to V{p), and 
the ehective masses in the plasma are dominated by vacuum and thermal masses (rather 
than coupling to the background p). Then the self-energies have the properties (33). In 


this case the behavior of 0 can be understood in the terminology of Brownian motion. 
Such systems have been studied by a number of authors [29H311EIHZ3], and it was found 
that the zero mode in this case follow a Langevin equation 


(46) 


PI + v|) .)>(«) + f ((- t')4>{t') = ^(«). 

Jo 


Here we have absorbed the local part of the r^-self-energy into a slowly varying on the left 
hand side, which therefore depends on T and p. We have also assumed that the dependence 
of t') out + t' is negligible. The validity of these approximations is restricted to the 
regime where the amplitude of 93 -oscillations is sufficiently small that the effective masses 
are not dominated by the 99 -dependent contributions. Otherwise the adiabaticity condition 


(27) is violated. This is consistent because we consider small oscillations near the minimum. 
Note that we do not need to require Tr^/br, -C 77 or 0^/0^ ^ H because the integral in 


(|46|) is cut off due to the hnite support of H^ {tp'). The size of this support is roughly 


given by the inverse of the (quasi)particle energies in the loop, but it is difficult to clearly 
identify the relevant scale in general m- The integral leads to dissipation, and ^s a 
Gaussian noise term with 


{pt)) = 0 and {i{t)pt')) = -n+(t - 


(47) 
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The relations (33) and (47) ensure that the fluctuation-dissipation theorem holds in (46). 


Strictly speaking the held in (46) is not the original held (j) on the contour C, but the 
combination 0c = + 0“)/2, where (j)'^ and 0“ denote the held with time argument on 

the “forward” and “backward” part of the contour, as indicated in figure [T] We ignore this 
subtlety here, as both notions give the same equation of motion for ip. 


The Langevin equation (46) can be solved analytically by use of a Laplace transforma¬ 


tion 


where 


0(t) = 0(t = O)A (f)-f0(t = O)A {t)+ A (t-T)^(0), 

Jo 


A (t) = / d^xA (x) 


(48) 

(49) 


can be analytically obtained by Fourier-transforming (34) in the Breit-Wigner approxima¬ 
tion as 


A"(f) ~ Z— 

a 6r) 


(50) 


By taking the expectation value (cf. ([^) of 0 in (48), we And (neglecting the po-dependence 
of Rell^, setting = Mn for the zero mode and neglecting terms suppressed by T^j/hlj,) 


p{t) = 


(p{t = 0)cOs{Mrit) + ip{t = 0) 




e ’"’'b 


(51) 


The solution (51) describes damped oscillations of p with the plasma frequency near 


the potential minimum. In spite of the integral in (46), non-Markovian effects do not play 


any role for p because T^ and Mn are in good approximation independent of p, and we have 
assumed that the temperature remains constant on time scales ~ l/T^. This is realistic 
if T, > H, and if the initial p is sufficiently small that the energy release into the plasma 
due to its dissipation is not significant. In this case an equation of the form ([^ holds, with 


m2 ^ 

V((p) ^ -fp\ 

r ~ r 

^ ip — ^ ri' 


(62) 

(53) 


That is, V{p) can be approximated by a parabola and is obtained from the free classical 
Lagrangian if one simply replaces the mass by the effective mass Mn- The damping 
rate of p is given by the thermal quasiparticle width of p; it is obtained by evaluating 
Imn^(p) at the quasiparticle pole po = Mn- Finally, each mode of 0 couples to the bath, 
but there are no interactions of the different modes with each other. This approximation 
only holds for small deviations from the ground state, i.e. small field values p. In cosmology 
this situation can be realized in the late phase of cosmic reheating after inflation. In this 
context, the temperature dependence of the rate F<^ has been studied in detail in eg. and 
it was found in eg that this temperature dependence can have a drastic effect on the 
thermal history of the universe during reheating. 
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For larger deviations from the ground state, the problem becomes much more involved 
because the quadratic approximation of the effective action near <^9 = 0 is not sufficient. If 
it is still possible to describe the dynamics in terms of an effective potential V(</9) and ([^, 
then V{(p) certainly must involve higher powers of if. This leads to “multiplicative noise 
terms” in the generalization of (46) [311 EIIj and the effective masses generally depend on 
f. Finally, we do not expect the simple relation (53) to hold, though it is frequently used 


in the literature without justihcation. In section 3.2 we show that one can nevertheless hnd 


an equation of motion of the form ([^ if changes slowly compared to the microphysical 
scales, but the expressions for V((p) and V^p are in general more complicated. 


3.2 Large field values: Slow rolling 

We continue to assume that all degrees of freedom except f are in equilibrium (or their 
thermodynamic state can at least effectively be characterized by an effective temperature), 
but we allow f to take large values, where the quadratic approximation ( |5^ of the effective 
potential is expected to break down. If 0 has only feeble interactions (which is required if 
the effective potential should be flat), then it is justihed to assume a separation of scales 


This implies that f changes adiabatically, and we can approximate 

f{t'Y ~ f{tY + n{t' - t)0(t)(^(t)”“b 


(54) 


(55) 


The change of the functions nc(t, t') and nc(t, t') with respect to the time coordinate t + f 


occurs on the same macroscopic time scales as f/f and H. Hence (54) ensures that (27) 
is fulhlled and we can expand nc(t,t') and in this coordinate analogously to the 


expansion (55). For the present purpose it is sufficient to keep the zeroth term in this 
expansion, i.e. to assume that there is no explicit dependence on t + F, meaning that 


licit, t') ~ nc(t — t') and licit, t') ~ nc(t — t'). This will be justihed later after (58). 


We can split the contour integral in the equation of motion (20) into 


n pcC)-\-ie 

dt' = / dt' + 

f C J—oo+ie 


dt'. 


(56) 


With (55), the only quantities under the integrals that depend on t' are nc(t — t') and 
nc(^ ~ t')- Since we are interested in physical times, t always lies on the “forward” part 


of the contour C. If t' also lies on the forward part (hrst term in (56)), then nc(t — t') is 


to be identihed with the usual time-ordered self-energy (often referred to as n”*"+(f — t') in 


thermal held theory). If t' lies on the backward part (second term in (56)), then nc(t — t') 
is to be identihed with a “Wightman type” self-energy without time ordering n^(t — t'), 
as t' is then always “later” along the contour than t, see e.g. [33 [62]. Their diherence can 
be identihed with the usual retarded self-energy, i.e. — t') — n”^(f — t') = n'^(f — t'). 


This allows to combine the two terms in (56) into a single integral 


dt'Ucit-t')= / dt'U^it-t') 


(67) 
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and likewise for nc(f — t'). The RHS of (57) is nothing but the Fourier transform of 
n^(f — t') evaluated at energy a; = 0. Some care should be taken at this point. The 


limit a; —)■ 0 is obtained due to the linear approximation (55), which is controlled by the 
separation (54) between macroscopic and microscopic time scales. It is probably best to 
think of cj as a parameter that characterizes this separation, i.e. u ~ where (p/(p 

should be identihed with the macroscopic scales H or T. With these considerations we 


hnally obtain a Markovian approximation to the equation of motion (20), 


[d'i + ml + {l- ud^) nj(-a;)] (p(t) + - + (1 - nj(-2a;)j (p(t) 


a ;=0 






= 0 . 


(58) 




Here we have given the retarded self-energies H;^ and 11^ an index to distinguish them 
for the self-energies for Xi. The parameter u should be thought of as much smaller than 
all microphysical scales and can be set to zero as long as a; < if < T^ < T, which always 
holds under the previous assumptions. Since the real part of the retarded self-energy is 
always symmetric and the imaginary part antisymmetric, we can easily identify 


lim (a;) 


a;^0 


lim Ren^(a;), 




lim 

oj—^-O 


R{ 


-UJ] = 


-i lim 


Imn«( 


OJ 


UJ 


(59) 

(60) 


and likewise for H^ 

With this limit and by comparing (58) with (1^, we can express the effective potential 
and dissipation coefficient in terms of H^and H^ The hrst line in (58) describes the evo¬ 
lution of (p in an effective potential with the quantum corrections to tne classical potential 
given by H^ and 11^: 




lim 


[ml + [1 


R{ 




-UJ 




A(/, -l- [1 


R( 




-2uj) 




— [i^l + n^(a;) |aj=o] T’ + ^ + n^(cu) |^=o 




(61) 


The second line describes dissipative effects, i.e. the damping of p due to the interactions 
with the bath of Ti-particles: 


T 


¥5 



(62) 


These results allow to justify why we could neglect their dependence on t -|- t'. The terms 
containing derivatives dt+t' acting on H^ and H^ are suppressed with respect to terms 


^From these expressions it is clear that all coefficients of y, and y in ( [5^ are real, 
effective equation of motion (58) is manifestly real. 


Therefore the 
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without derivatives due to the separation of macroscopic and microscopic scales. For 
instance, is suppressed by three small parameters, two powers of a coupling 

constant and a time derivative. All terms in the hrst line of 


which give the effective 
potential (61), contain two or less small parameters, hence can be regarded 

as a higher order correction. In spite of this, we have to keep the terms in the second 


line in (58), which in principle are also suppressed by a time derivative and two powers 
of a coupling constant. The reason is that these form the leading order contribution to 


the dissipation coefficient (62), while in the effective potential (61) they would just act as 


higher order correction. In this sense, (58) is a consistent expansion in gradients and the 
small coupling constants. 


These results (61) and (62) should be compared to (52) and (53). One obvious (and 
expected) difference is the appearance of higher powers of (p in V((^). Another crucial 


difference is that the dissipation coefficient (62) is evaluated at a; = 0 and not ai u = M, 


‘■ni 


as in (53). This supports the interpretation of u: as the rate at which (p changes: During 


oscillations around the minimum, this rate is given by the plasma frequency Mr^] if the 
held slowly rolls down a hat potential, then this rate is much smaller than all other scales 
(practically a; ~ 0 within loop diagrams). 

Indeed, if we drop the terms oc in (20) and use (56) and (57), we hnd 


( 5 / + mlMt) + / = 0 . 


(48), 


This equation is formally the same as (46) with ^ = 0, so we can read oh the solution from 


Lp{t) = = 0)A (t) + ip{t = 0)A (t). 


(63) 


Again using the Breit-Wigner approximation for (34), we can bring this into the explicit 


form (51). Hence, we have recovered the behavior of Brownian motion for small ip from 


the more general expression (20). 


4 A simple scalar model 

In the following we illustrate our results in a simple scalar model, in which the “bath” 
consists only of one other real scalar held x. 


C = p„4>8^4> - -y4> - X - - ^x - X ■ 


(64) 


In order to obtain explicit expressions for V{p) and T^ from (61) and (62), we have to 
evaluate the loop diagrams given in hgure and 


Here we focus on the case of adiabatically changing large p discussed in section |3.2 
For the case of oscillations near the potential minimum, F^ has been calculated in detail 
in 
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Figure 2: Self-energy diagrams contributing to mass correction at leading order. Solid line 
corresponds to p and dashed line to x- Note that the functionals 11^ and 11^ should strictly be 
represented by amputated diagrams. The external lines are added for representative purposes 
and in the effective action they are to be replaced by (p. 



Figure 3: Feynman diagrams contributing to 11^. Solid line corresponds to rj and dashed line 
to X- 
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4.1 The effective potential V{if) 


From (59) it is clear that V{if) is obtained from the real part of retarded self-energies. 


To establish a consistent perturbation theory, these have to be evaluated with resummed 
thermal propagators. For the present purpose it is sufficient to evaluate these in the pole 


approximation of (39) for r] and y. The potential in (64) in terms of r] and x reads 

{ml + hif^/2) 


y{v,x) = 


m. 


Aw 




with ml = ml+ 




4 I 2 I "'„2, 2 


h 


^-vx +jVX 


(65) 


4.1.1 The self-energy 11^ 

The real part of n^|^=o is dominated by the well-known local diagrams due to the vertices 
and hrfx ^see figure]^ It reads 


= 


Xd 


Xd 


1 r „ M h 




[l + 2/B(fiJ] 




Ml 


{AtiY 

f dp^ 


[iV. + T(M,^/^^)] + 2(4,)2 


(27r)3 2X1^ 

[N, + T{MI,pY] 


X(j) 


with 


) (27r)3 

T(M2,p2) = 1 




2M 


2 

log^, N^ = - + loglvr 
jjL^ e 


7, 


( 66 ) 


(67) 


where we have used dimensional regularization. The divergent part {N^ terms) can be 
eliminated W renormalisation of and coupling constants at zero temperature, using 
MS schemejj In the following we assume that this has been done and that vacuum masses 
and coupling constants are the renormalized parameters at T = 0. In this work we are 
mainly interested in thermal corrections, which are most relevant if the temperature is 
larger than the effective masses (otherwise the thermal corrections are negligible since 
they are Boltzmann suppressed due to the appearance of the distribution functions fs in 
the thermal propagators). For ^ T, the finite part (i.e. the T(M^,/i^) terms), 

remaining after subtraction of the divergences, can be ignored since it is subdominant 
compared to the /s-part (the third line of ([6^) Due to the momentum-in dependence 

®Note that the temperature-dependent divergent terms, which come from the thermal corrections con- 

are sufficiently taken care of by the temperature-independent 


tained in and multiplying N, 


counterterms from and h in the third line of ( 661 , which in fact represents the thermal corrections. 

^This can be seen as follows. Let’s choose /r to be /r = Mi = Max[Mjj, evaluated at initial 
value of y. Then < Mi for all the times. Now we will show that jA^ T(M^, )| <C \rj)T^ 

and hT"^ when M^,M^ T (so Mi T). To this end it suffices to show 

A/T‘2 

This is clearly the case since, for example, 
1 (note that |eloge| < e~^ for 0 < £ < 1). 


j\yf2 

X^ Ml log 

^ A^T^ and 

A/f^ 

h Ml log 

< 


log $ = 

log^ 
Mf Aff 

< Mf 

for 

< < 
Mf ^ 
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there is no need for a finite T wave fnnction renormalisation. For <C T, one can 

perform the integrals analytically, for example, 


f dp^ h 

J (27r)3 2f]^ 


/n(^x) 


rs> 


A 


rji2 


The same applies to the diagram with r^-loop. We thus obtain for -C T 


( 68 ) 


„R (/i + A^) ; 

«■ - 24 


(69) 


There are also contributions from the setting sun diagrams in figure but they are of 
higher order in h and X^. 


4.1.2 The vertex functional 11^ 

We now turn to the vertex functional, which is given by the “fish diagrams” a) in figure 
1^ The /^-independent part of these diagrams is again divergent; for vanishing external 
momentum it is given by (including the analogous fish diagram with r^-loop) 


2(47r)2 


W - log 


All 

/r2 _ 


+ 


A? 


2(47r) 


W - log 


/r2 _ 


(70) 


Since we are only interested in the case of vanishing external momentum, we can simply 
absorb the W terms into the renormalisation of the vacuum couplings h and A^. The 
remaining hnite part (log terms) can be ignored for T in comparison with 

/s-part given below, in the similar way to the case of 11^,^ 

The /s-part from the fish diagrams is given by the the standard expression (see e.g. 
section 4.2.2 in ra). Let us focus on the X"loop with T S> M^, 


TT^I — 

-*■-*■05 |aj=0 


47r 


dVL, 








/b(A) “ 


A^ 47r 
4 (27r)= 


dVL, 


x/ai, - 


Mn 


V 


fsiXlr,)- 


(71) 


The above integral can not be done analytically, but it is easy to see that it is strongly 
dominated by small 12 when T. Therefore we can approximate /b(12) ~ T/12 

by expanding the exponential to get 


h2 dvr 
T (27r): 


[ dn^^ 

'Vfv 




M2 

—/n(^^x) ^ 


h2 dvr 


T./U 

da ^ 


Ml 


‘'X 


'Mv 




h2 T 
327r Mv 


(72) 


^Setting again /x = Mi = Max[M^,M^] evaluated at initial value of ip, it is straightforward to see that 


•1 

IoSmT 

be shown for M~ 


iln Ihi Ino- lua 


< ^ <C where we used the same tricks as in footnote The 


same can 


X- 
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The contribution from the r^-loop can be obtained from the above expression by replacing 
X V and h —)■ and altogether we obtain for T 


nJU=o ^ 


T 


\2 


T 


327r Mv 327r M„ 


(73) 


4.1.3 The full potential 

If the temperature is much smaller than the effective masses, then thermal corrections are 
negligible. In this case radiative corrections are mainly interesting if some mass squares are 
negative, leading to symmetry breaking; otherwise they just lead to small modifications of 
the potential. If T is comparable to the effective masses, then the loop integral in general 
cannot be solved analytically. For the interesting case M^, <C T, we can use the above 
results in (61) and hnd 


= (^ml + 


[h + A(^) 


24 


rp2 


If. T XI T 

+ 777 ( ~ 


3! 


327r Mv 327r Mr, 




(74) 


In order to check for which range of field values the assumption M,j,M^ T can be 
justified, we can estimate the effective masses as 


~ + (^ + 


+ (^ + ^x) T7 + 77V’ 


T2 


^2 

2 
h 


(75) 

(76) 


This leads to the conditions 


< 




Arf 


, 79 " < 


1 - 


(fe+Ay) 


24 


T^-ml 


h 


(77) 


These conditions imply that the (p-dependence of the effective masses in (74) is not neces¬ 
sarily negligible. Taking this into account, we integrate (74) and obtain 

T 


V(v) = 

where 




+ 


4! 


+ 


2887r 


4(M;; + M^^ - + hMj , (78) 




Mr,{ip = 0)^ 

= 0)^ 


T‘2 

+ (^ + A^) 

rj-^2 

+ (/^ + 


(79) 

(80) 


Figure]^ shows the effective potential (78) as a function of 79 for various temperatures. 
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v(ip)lm^‘ 




Figure 4: The full effective potential V(ip) as a function of tp for various temperatures, normal¬ 
ized to (figure a)) and to zero temperature value Vo(<p) = V{ip)\t=o (figure b)). Here we 
set X,p = = 1/5 and h = 1/30. Different values of temperature are indicated 

by the color coding with the pure blue line being for the lowest temperature T = 2Om0. When 
the color varies from blue towards red, the temperature increases with a step AT = 20to^. 
The dotted line is used to denote the regime where the condition < T is not satisfied. 


4.2 The damping coefficient 


Damping coefficient (62) depends on ip through the effective masses and has to be evaluated 
at a; = 0. 


4.2.1 Contribution from 


Using (60), we can determine the coefficient lim(^_,.o by calculating Imn^(p) from 




diagrams given in hgure for vanishing p. The leading order contribution comes from the 
setting sun diagram b) in figure since diagram a) gives vanishing result for kinematic 
reasons]^ The calculation of the loop integral is rather technical and is summarized in 
appendix [A} For T, the result reads (cf. (122)) 


dM 


Ri 




, Vm, 


(81) 


4.2.2 Contribution from the vertex 11^ 

Let us focus on diagram a) in figure which gives dominant contributionp^ The contri¬ 
bution from the analogous diagram with r^-loop can be obtained by proper replacements 


® There are two kinds of processes contained in cuts through the diagram in figure a) for F,^ = 0, 
decays cj) —>■ b/’/’ scatterings ^ (and their inverse). The former are clearly kinematically 
forbidden. For vanishing external four-momentum, the latter effectively correspond to 4>4> decays 

and inverse decays, which are also kinematically forbidden on-shell. Finite width corrections are of higher 
order in this case. 

for the argument that the vertex diagram b) in figure is subdominant. 


^°See section 


4.2.4 
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b) 




Figure 5: Feynman diagrams contributing to 11^ with only ^-propagators in loops. Dashed 
line corresponds to x- Note that there are analogous diagrams which involve rj-propagators 
in loops. 


In the zero-width approximation (39), the imaginary 
^ and diverges for to = 0,^ One can 


of conpling constants and masses, 
part Imll^ vanishes for any hnite positive to < 2M- 
nnderstand this from kinematic reasons when recalling that the imaginary part can be 
interpreted in terms of microphysical processes when applying the optical theorem at hnite 
T IZg , see hgnrej^ A non-zero can be obtained by inclnding the hnite width of 

the x-propagators in the loop, which parameterizes the ehect of scattering processes with 
more vertices [631 [75] , see hgnre This aspect has been pointed ont in [ZZIEHI. 

The loop in the “hsh diagram” a) in hgnre can be expressed as 


ImIlJ((o) 




2 J (27r 


^ ^ ^iA<(po)A<(co-po) 


X 

. <1 


r d^p 

] W. 


4 (1 + /b(Po) + /^(a; - po)) Px{Po)pxi^ “ Po)- 


(82) 


For notational simplicity we continne to snppress the dependence on spatial momentnm; 
the external spatial momentnm is vanishing, and for all qnantities nnder the integral it is 
given by the loop momentnm. Here we have hrst nsed (40) and the KMS relation (33) 
to express ImH^ in terms of H^. Using the Feynman rnles sketched after eqnation (45), 
which are derived in detail in 1371. this fnnction can be calcnlated from an integral over 


alone. With (|41|), this can hnally be expressed as an integral over a prodnct of spectral 


densities, which is particularly easy to be evaluated in the approximation (39). In the 


analogous expressions (94) and (123) for the setting sun diagrams one can indeed use the 


approximation (39). In (82), on the other hand, we cannot neglect the hnite width and 


These two facts can be seen from figure [l0|b) and c) , whic h depict the supports of the spectral densities 


in the expresion (82) for Imll^, as discussed in section 


5.2.2 
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have to use the pole approximation to 
Px(Po) - 


, see e.g. ua for a discussion, 

-2Imn^(po) 


(Po - ^x)(Po + ^x)(Po - ^t)(po + ^t) ’ 


(83) 


where = f2^ + |r^. For a; <C F^ the quasiparticle resonances at po — ^x Po — 


of (83) in (82) overlap in the entire integration region, and the integral is strongly dominated 


by this “pole contribution”. Inserting (|83|), we obtain 

d\v 




2(27r) 


-p 2Re 


Imn^(a; - (2^) (1 + /i?(^x) + /b(^ “ ^x)) 


La 






Imng(a; + ~ + ^x)) 

(a; + Q*y - Ql) + (lmn^(a; + Q*) 


(84) 


We now take the limit a; —)■ 0: We insert (84) into the relation lim^_,.o 9tjli^(a;) = 


lim^_,.o with approximations 


Imn^(a; — f2^) ~ Imn^(a; — f2^) , Imn^(a; + ) ~ Imn^(a; + f2^) 

and then replace Imn^(a; ± f2^)|tj=o —t in the denominator to obtain 


(85) 


lim 5^n^(ci;) = lim 


ilmn^(a;) 




cj^O ijJ 

40^2 f 




—i 


tt^T 

40/^2 


d\p\ 

d|p| 




(F4 + 68 F 2 f 22 + 256(2^ F^ (cosh(f2 JT) - 1) 
P^f22 


256f2^F^ (cosh(f2^/T) — 1) 


( 86 ) 


The calculation of the momentum-dependent thermal width F^ from the two loop diagram 
in £gure[^is generally difficult. The full F^ is sum of contributions from diagram a) and 
b) in hgure 1^ i.e. F,^ = F^“^ -|- F^^\ The calculation of the imaginary part of diagram a) 


in hgure 


(cf. (164)) 


6 is presented in detail in appendix 


B 


and Fx“^ can be approximated <C T as 


(a) T’’^ 

F(“) ~ ^ with 7 ^“) = 

^ 2567r3 


(87) 


The main momentum dependence is due to the time dilatation; the factor is due to 
the extended lifetime of a relativistic X"Particle in the bath rest frame. The contribution 
from the diagram b) in figure should have essentially the same behavior as that of a) as 
long as <C T because the loop is dominated by hard momenta ^ in this 
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Figure 6: Feynman diagrams for self-energy of x Spying finite width of the y-P^opagators in 
the loop in figure [^aj. Dashed line corresponds to y and solid line to rj. 


case. Therefore one can obtain the contribution T^^ by replacement —)■ h in (87) and 

by taking into account a proper symmetry factor of the loop: 


~ ^ with 7 ® = 


(b) 


n 


Sh^T'^ 

2567r3' 


( 88 ) 


The integral then reads 


lim 9^n^(a;) -i 


ui^O 


AOnh^ 


T3(A2+3h2) 


d|p| 




(cosh(r2^/T) — 1) 


(89) 


This integral is strongly dominated by the region T. We expand the integrand in 

r2^/r up to second order and integrate from up to = T to obtain 




807rh2 


| a ;=0 


—I 


log( ^ 


T(A2 + 3h2) ^ \M, 


(91) 


This is for the fish diagram a) in figure]^ with y-loop. The contribution from the analogous 
hsh diagram with 7 -loop can be obtained from above expression by the replacement —>■ 
and / (A^ 3h^) -)• A|/(A^ 3h^). 


should be pointed out that for 0 < w < F^^ and with the previous assumptions, the integrand in 
(841 formally exhibits a sharp peak at = r2/(8a;). In this region the integrand can be 


approximated by the Breit-Wigner curve. In the zero width limit the contribution from this region can be 
estimated as 


(2 + 




v2/3 


(7^a;)2/3 cj4/3 


) 


30727rT (cosh 


- 1 


(90) 


Though this contribution smoothly vanishes in the limit a; —>■ 0, it is much bigger than (911 for most finite 


0 < w < F^ even if w F^^. This appears to be at odds with our intuitive argument that the result of 
(841 should be independent of w if |a;| <C F^^. The contribution (90) is not physical and we artificially 


introduced by the approximation (85), which is only valid in the region 12 
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4.2.3 The full damping coefficient 


Using above results in (62), we obtain the full damping coefficient for T 


T ^ 


(47r)3M^ 




SOvr 




Xl + 3h2 

X. 


1 + 


A2 


A2+3h2 




which is plotted in figure as a function of p for various temperatures. 


( 92 ) 


4.2.4 Physical interpretation of the coefficients 


The parametric dependence of (81) can easily be understood physically; the rate for 2 —2 


scatterings is proportional to and grows oc due to the higher density of scattering 


partners, just as a quasiparticle damping rate would. The coefficient (91) shows a more 


unusual behavior: most notably it decreases with temperature and with A^, both of which 
may seem counter-intuitive if one has the scattering interpretation in mind. It should, 
however, be pointed out that neither the limit T —)■ 0 nor A^ —)■ 0 can be applied to the 


approximate analytic formula (91) because it is only valid under the assumptions ikU < T 


and (p/p < oc A^, in which case the limit a; —)■ 0 is only justified. In addition, the 
limit A^ —)■ 0, which makes oc A^ vanish, is not allowed for following reasons. In 

the approximation F^ = 0, the cut through the diagram in figure at finite T includes 
various processes, which are obviously kinematically impossible for any finite uj < 2M^ 
as explained in the caption. Moreover, for a; = 0 the rate has an unphysical divergence, 
which can be traced back to the infrared divergence of the Bose-Einstein distribution for 
the quasiparticlesj^ The divergence in the limit w —)■ 0 is regularized by the finite width 
F^ 7 ^ 0 of the intermediate y-particle in the contributing processes, see figure Therefore 
F^ needs to be kept non-vanishing. 

Though not infinite, the cross section is resonantly enhanced with finite F^^ as l/F^^^, 
and it becomes larger for smaller F^ oc A^. This explains the A^ in the denominator of 


(91). The resonant enhancement of the “self-energy correction” in figure]^ a) is not present 
in the “vertex correction” in figure |^b). This allows us to neglect the vertex correction, 
which is of the same loop order as the self-energy correction. Quantitatively there is not 
much difference between the rate at finite a; F^ and w = 0, as one can see from figure [TOj 
The resonant enhancement would be weaker if we were dealing with fermions with gauge 
interactions in the loop (instead of y). The reason is that the damping F^ in the scalar 
model only appears at two-loop order (i.e. from the setting sun diagram in figure]^, while 

we were dealing with real quasiparticles in the initial and hnal states (rather than those with 
vanishing external four-momentum in the condensate), then this divergence would be regularized by the 
thermal mass; there would be no such case as w = 0 because the self-energies are evaluated at the 
quasiparticle mass shell. In a plasma the quasiparticles at rest are massive due to the thermal mass even if 
their vacuum mass is zero, such as for photons. In the limit T —>■ 0 the thermal mass disappears and only 
the vacuum mass remains (which can be chosen to be zero), but the divergence from the Bose-Einstein 
distribution also disappears. 
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b) 






Figure 7: The damping coefficients as a function of ip for various temperatures, normalized 
to m,f,. In these plots the choice of the parameter values and the color coding of lines for 
different temperatures are the same as in figure [^except that the dotted lines here denote the 
regime where the condition < T is not satisfied. Figure a) shows the full damping 

coefficient F^. In figure b) we plot the ratio ri/r 2 with Fi and r 2 being contributions from 
(setting-sun diagram) and 11^ (fish diagram), respectively (i.e. the first and the second 
term in (92)). This plot shows that Fi is dominant only for very small p and otherwise r 2 is 
the main contribution to the damping. Fi and F 2 are separately plotted in figure c) and d). 


27 









the mass correction is of order (from the tadpole in figure]^. This makes T^/hl^ very 
small for the scalar model. 

Finally, in section II of [72] it has been claimed that all dissipative effects vanish if 
the light held y is in a vacuum state. This appears to be conhrmed by our result for fl^. 
Physically this conclusion seems, however, counter-intuitive. The basic laws of statistical 
mechanics imply that the energy in an interacting physical system should relax to a state 
of equilibrium, which implies equipartition of the energy amongst all degrees of freedom. 
If dissipation were absent for T = 0, then all energy would forever remain trapped in (p if 
the initial state of the system contains no y-particles. The reason why we hnd vanishing 
dissipation rates for T —0 in the present calculation lies in the approximation a; = 0. 
Note that Imn^(ci;) in (84) gives a non-zero contribution for hnite uj and T = 0. The 
precise determination of the physical damping rate at T = 0 would, however, require some 
extra work and a more consistent treatment of renormalization. 


5 Discussion 


5.1 Main results 

Effective equation of motion for the field ip - In section we have shown from hrst 
principles that the expectation value (p of a scalar field in a medium approximately follows 
the Markovian equation of motion (58), which is of the type ([^, if the coupling to the 
medium is weak and the properties of the medium change sufficiently slowly. That is, the 
system can be characterized by a complex valued effective potential. Its real part V((p) is 
the usual finite temperature effective potential, and the imaginary part is the dissipation 
coefficient. They depends on temperature and the held expectation value (p. The equation 


of motion (58) and the general expressions for V{ip) in (61) and T^ in (62) are amongst 

^ isi, 


The expression (62) extends our calculation of in 


the main results of this work, 
which is valid near the potential minimum, to the case of large held values in the slow-roll 
phase. They appear to be consistent with the equation of motion for the held (j) found in 
the literature (see e.g. [HU |3l| and references therein). However, our derivation is much 
shorter, as we directly seek an equation of motion for ip. Together with the expressions 
(28) and (29) for the propagators, the equation of motion (|^ with V((p) in (61) and T 


m 


(62) allows to describe the nonequilibrium dynamics of a slow-rolling scalar held entirely 


in terms of Markovian equations. 

For small values of (p, and if the medium is composed of a sufficiently large thermal 
bath, it is well-known that the held 0 is exposed to Brownian motion, while its expectation 
value ip performs damped oscillations and relaxes to its minimum on a time scale that is 
given by the thermal quasiparticle width in the plasma. We have recovered this behavior 
as a limiting case in (63). Far away from the potential minimum the behavior is very 
diherent. 
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Figure 8 : Cut through the hsh diagram with zero width for the propagators. Solid lines 
correspond to (j) and dashed lines to x- Due to the optical theorem at finite T, the imaginary 
part of the hsh diagram can be interpreted in terms of physical processes in which the external 
legs and cut propagators appear as initial and hnal states The diagrams on the right 

show those processes in which initial and hnal state both contain particles, as only these can 
contribute in the present calculation. Note that all y-purticles here are on shell. Diagram a) 
corresponds to the process (jicj) —>■ xXt which two (j)-quanta out of ip-condensate become two 
X-quanta. If the total four-momenta of the two (p-quanta is (w, 0 ) with co < 2M^, this process 
is kinematically not allowed for on-shell x-partides. Diagram b) represents the scattering 
x4> xi' between a x-quantum and a single p-quantum in the condensate (in both, initial 
and hnal state). If the sum of four-momenta of the p-quanta is (u, 0) with any w 7 ^ 0, then it 
is clear that such process is also forbidden for on-shell x-ps-rticle due to the energy-momentum 
conservation. Diagram c) represents a decay process x x4"P- Here the initial state contains 
a single x-Quasiparticle with non-vanishing four-momentum. In the hnal state, there is also a 
X-quasiparticle and two (p-quanta with combined four-momentum (uj, 0) have been added to 
the (fi-condensate. This process is kinematically forbidden on-shell for any to 0, but the cross 
section is divergent for w = 0. Similarly, in diagram d) an p-quantum in the (p-condensate 
“decays” into an p and two x-Quasiparticles. The resulting total cross section is vanishing for 
any hnite uj < 2M^ and diverges at lo = 0. These intuitive interpretations hold at the level of 
individual quanta: The (p-condensate is a superposition of inhnitely many p-quantum states 
with different particle numbers. However, if one decomposes the condensate into eigenstates 
of the free Hamiltonian (with well-dehned particle number), then the same argument can be 
drawn for any multi-particle state. 
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Figure 9: Cut through the hsh diagram with a Bnite width (T^ ^ 0) for the x-propagator 
in the loop. If one evaluates this hsh diagram, then it includes contribution from process 
VXX VXXy which is mediated by an off-shell intermediate x-quasipartide. This process is 
allowed and hnite for w = 0. 




Figure 10: Support of the spectral densities Px(Po) stnd PxiPo — in po-|p| plane. The solid 
blue and red lines correspond to the mass-shell Pq — |pp = and (po — uj)^ — |pp = M^, 
respectively. The dotted lines represent the approximate support of the spectral densities if 
the hnite widths F^ are included. Figure a), b) and c) correspond to the cases ofoj> 2M^ 3> 
F^, 2Mx > uj > Tx and 2Mx ^ F,^ > w, respectively. In figure a), the region where two 
mass-shell curves intersect gives rise to the non-vanishing contribution from kinematically 
allowed processes of the on-shell x-particfes, such as decay and scattering. 




Figure 11: Spectral densities Px(Po) (blue line) and Px(po — (^sd line) as a function of po 
for hxed |p| (e.g. for |p| = OJ. Figure a), b) and c) represent the same cases as in figure [I^ 
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Effective potential and dissipation coefficient in a scalar theory - In section 
we have applied the formulae derived in section 3^ to a simple scalar model. Analogue 
results for the small held case discussed in section 3.1 are given in [75]. The hnite temper¬ 


ature effective potential and damping coefficient are approximately given by the analytic 
expressions (78) (equivalently (74)) and (92). Our results appear to be in some tension 


with the previous literature, which we discuss in section 5.2 


Loop integrals at finite temperature - In the appendices we have provided general 
expressions and analytic estimates for nontrivial two-loop integrals for setting-sun diagrams 
at hnite temperature. This includes the setting-sun diagram at vanishing external four- 
momentum (appendix I A. 2 D and thermally on-shell cases (appendix |B.2[ ). In appendix]^ we 
have given analytic results for the angle integrals in setting-sun diagrams. These formulae 
will be very useful for further calculations in hnite-temperature scalar held theory. 


5.2 Comparison with previous results 

Our results for the damping coefficient T^ diher from those obtained in the previous liter¬ 
ature in two ways. 


5.2.1 cu = 0 versus u = Mr, 


A main diherence between the damping coefficient (53) during oscillations near the poten¬ 


tial minimum and the expression (62) for large held values in a slow-roll phase is that the 


diagrams in the former are evaluated on the quasiparticle mass shell, while in the latter 
they have to be evaluated at vanishing external energy (i.e. at vanishing external four- 


momentum, which is oh shell). In spite of this, the coefficient (53) is frequently used in the 
literature in all regimes, which is clearly incorrect. This point has previously been realized 
by some authors, see e.g. [HU [771 Elj- In spite of this realization, the authors of [211IM] 
used the on-shell damping rate (53) in their calculations, probably for simplicity. The 
authors of mi, on the other hand, calculated with vanishing external four-momentum. 


The fact that the damping coefficient in (62) is evaluated at vanishing external four- 


momentum raises the questions to which degree our results are sensitive to various known 
infrared problems of thermal held theory [551 - 157] . and whether further resummations and 
the inclusion of vertex (and ladder) diagrams are necessary to obtain physically consistent 
results. For example, in the calculation of 11^ we neglected the vertex diagram shown 
in hgure [^b). This can be justihed in the present case because it is not subject to the 
resonant enhancement l/T^ of the hsh diagram in hgure a), cf. (86). In the calculation 


of 11;^, the main argument in favor of our approach is that the thermal masses in (75) - (76) 


already appear at order and \/h, while the thermal width is of higher order 

in the couplings. This tends to suppress all sorts of hnite width and vertex corrections. 
However, this handwaving argument is of course not a strict proof, and we postpone the 
clarihcation of these issues to future work. 
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5.2.2 Evaluation of the fish diagram 


The dissipation coefficient at leading order is given by the imaginary part of the retarded 
self-energy from the fish diagram, Imn^(a;). The fact that it should be evaluated at 
vanishing external four-momenta (i.e. a; —)■ 0 in (82)) was previously pointed out in 
[HlEH]. However, the results obtained there at first sight appear to differ from ours. In 
the following we try to understand the origin of the difference. Again recalling the optical 


theorem, we can identify the region in the integration volume of (82) where both spectral 
densities are thermally on-shell (i.e. at the intersection of two mass-shell curves in figure 
10 a)) as the non-vanishing contribution from kinematically allowed processes (e.g. decay 


and scattering) of the on-shell y-quasiparticles. From the kinematic considerations given 
in figure it is clear that such a region exists for uj > 2M^. For F^ <C 12^ it strongly 
dominates the integral (82), see e.g. |63l [75] for a detailed discussion of this integral. In 
this case one can use the zero-width approximation PxiPo) — 27rsign(po)(5(po ~ 
evaluate the integral 


For F^ < a; < 2M^ (see figure 10 b)), the use of this approximation 


is not allowed. In this case the integral is dominated by the regions in which one of the 
Px is on-shell (see figure 11 b)). Within these regions, the part, where the distribution 
functions Jb have their maxima, gives the biggest contribution. For the case w F^ 
under consideration in this work, the pole regions overlap in the entire integration volume. 


as depicted in figures 10 c) and figure 11 c). This leads to a strong enhancement, and the 


integral is generally dominated by this pole region po — ^x- 


expression (83) for the spectral densities to obtain the results (84) and (86). 


The setup discussed in [77117SI is an interesting exception, in which the integral ([M|) 


for cj Fy is dominated by the off-shell regions in spite of the alignment of the two 


quasiparticle peaks illustrated in figures 10 c) and figure 11 


contribution from the poles becomes Boltzmann suppressed for ^ T. 


The reason is that the 
For sufficiently 

large M^, the integral in (82) for a; —)■ 0 is dominated by the po — 0 region. In the model 
under consideration here, all thermal corrections to the y-propagators become suppressed 
in this limit, and one essentially recovers the vacuum limit. By adding another fields a 
with effective mass Mo- T ^ M^ and an interaction term the authors of [ZZIEH] 
have constructed an interesting scenario in which the loop integral is dominated by off-shell 
regions, but thermal corrections are not negligible. If the y-self-energy is dominated by 
the interactions with a, then its temperature dependence is crucial and affects F,^. This is 
the origin of the qualitatively different behaviour of the results found in those articles, as 
compared to ours. The different regimes are discussed in more detail in [78] . They appear 
to be consistent with our results where the range of applicability overlaps. Practically the 


integral (82) can be evaluated using the approximation 


Px 


2Imn(^ 


(93) 


This approach was referred to as the low-momentum approximation. Unfortunately this 
introduces a strong sensitivity to the infrared behaviour of the self-energy, which for a ycx^ 
coupling is rather complicated and not fully understood [38]. The strongly hierarchical 
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arrangement T and T, corresponds to a rather specific corner in parameter 

space, but is cosmologically very interesting because it has been argued [HE] that it provides 
a viable “two-stage mechanism” for warm inflation [55] , 

6 Conclusions and outlook 

We have studied the effective action and equation of motion for the expectation value 
93 of a scalar field 0 in a dense medium from first principles of nonequilibrium quantum 
field theory. We focused on two cosmologically important processes, damped oscillations 
near the ground state and a slow-roll phase in a fiat potential. In a series of controlled 
approximations, we showed that these processes can be described in terms of Markovian 
effective equations of motion for both, the occupation numbers in the plasma and the field 
expectation value. This allows to describe the system in terms of a complex valued effective 
potential. The real part V((p) is the usual effective potential at finite temperature, which 
includes radiative and thermodynamic corrections to the scalar potential. The imaginary 
part is the dissipation coefficient that leads to damping and particle production. 

When applying our results to the damped oscillations near the potential minimum, 
we found that our method reproduces the well-known results that have been obtained 
by studying the Brownian motion of 0 near the ground state. As expected, the effective 
potential and damping coefficient for 93 coincide with the thermal mass and width of quasi¬ 
particles in the plasma. However, far away from the minimum we found a very different 
behavior. In particular, during a slow-roll phase, the loop integrals that determine the co¬ 
efficients in the effective action have to be evaluated at vanishing external four-momentum 
and in resummed perturbation theory. This shows that the common practice, which uses 
the thermal mass and width of quasiparticles as the order-of-magnitude estimates for the 
thermal corrections to the effective potential and dissipation coefficient, is clearly incorrect 
in this situation. 

We illustrated our results in a simple scalar model. In this context, we provided explicit 
expressions and analytic estimates for the setting-sun and fish diagrams at finite temper¬ 
ature. These will be very useful for future computations in scalar field theories at finite 
temperature. 

Our results mark an important step forward towards a quantitative understanding of 
the evolution of scalar fields in the early universe. For instance, they can be applied to study 
the fate of moduli, curvatons, axions and other scalars with a fiat potential during and after 
reheating. They may also be used to models of warm inflation. Finally, the scalar field 
may also represent an order parameter in applications outside the domain of cosmology. 
Our analysis is, to the best of our knowledge, the most systematic and comprehensive one 
of the subject to date. However, it still relies on several restrictive assumptions. We have 
used thermal propagators in Minkowski space to evaluate loop integrals. This can only be 
justified if the particles in the loop have reached a local kinetic equilibrium, and when the 
thermal masses in the plasma are larger than the rate of Hubble expansion. Both of these 
assumptions should be relaxed in future work. Finally, the model in which we illustrate our 
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Figure 12: Setting sun diagram for interaction. Solid line is rj, while dashed line being 
X- 

results is a simple toy model in which the primordial plasma is composed only of another 
real scalar held. In most realistic applications, the scalar held in question couples directly 
or indirectly to fermions with non-Abelian gauge interactions. Some progress towards 
an inclusion of medium ehects on the ehective potential has recently been made e.g. in 
[Ml ESI El ESI [m [721 ESI Enra , and phenomenological implications have been discussed 
in [751 IMbl02j , but it is still a long way to go to a complete quantitative understanding of 
scalar helds in the early universe in realistic models. 
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A The self-energy from setting-sun diagram 


The main goal of this appendix is to evaluate the self-energy in (|60j) from the setting-sun 
diagram shown in hgure 12, using the zero-width spectral density (39) neglecting 
For earlier discussions of these diagrams, see e.g. [TUmilT] . We consider the spectral 
self-energy H”, which is related to the imaginary part of the retarded self-energy. 
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A.l General expression 

The general expression of n~p(po) = 2*Imn^p(po) is given by 


n^p(po) = —^ 


ih'^ f d^q d^k dH 


(27r)4 (27r)4 (27r)^ 


{2Tr)^6^^\p -q-k- l)p^{q)p^{k)p^{l) 


X 


(1 +/b(%)) (1 +/s(^o)) (1 +/s(^o)) — /b(?o)/s(^o)/s(^o) , (94) 


which can be derived analogously to (82). After performing an integral over q for p- 
propagator, we obtain (suppressing subscript p for If”) 


np(po) = -^(T>p(po) + Spi{po) + Sp2{po) - Vp{-po) - Spi{-po) - Sp2{-po)) (95) 


with 


^^p(Po) 


‘5pi(po) 


‘5p2(Po) 


Here 


TT .J {Po) / dill / df2kArf(f2k,fii,-Ai)X(ai,6,c), 

8(27r)5 




TT 


8(27r) 


dkli / dhlk-^si(f2k, f2i, Ai) X(ai, 6, c), 


'Mv 


27r 


8(27r)5 


S{Po) 




-^s2(^k7 ^1? -^ 2 ) ^(^2; c) . 


(96) 

(97) 

(98) 


Fd{flk, hli, A) 
Xsi(r2k, fli, A) 
-^s2(f2k, f2i, A) 
ilf 
kl2 

Ai 

A 2 

X(a, b, c) 
z 

I{x,y,z) 


(1 + /^(f^k)) (1 + /^(f^i)) (1 + /s(A)) — fsi^k) f b{^\) f b{.A), (99) 

(1 + /B(a))(l + fBm)fB{A) - /B(f2k)/B(f2i)(l + /b(A))(100) 

/B(hlk)(l + /^(hll)) (1 + /b(A)) — (1 + /B(f 2 k))/B(f 2 i)/B(A)( 101 ) 

max [po — rii ± Mj , (102) 

max [rii — Po + Mrj, , (103) 

hlk + 111 — Po, (104) 


flk 111 T Po, 



d{I{x,y,z)) 
\/l{x,y,z) ’ 


a + by cx^ 


{l-x^){l-y^)-{z-xy)\ 

yl2_ 2_k2_l2_M2 

— -rrrTTi- - j 


1,2, b 



(105) 

(106) 

(107) 

(108) 

(109) 


The V terms represent decay and inverse decay p ■H- pxX- Note that the integration limits 
in (96) imply that Xp(po) vanishes for po < 2M^ + Mri, as expected from energy-momentum 
conservation in the processes p ■H- pxx fh^'f cire impossible for on-shell particles. 8Spi(po) 
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and 5p2(po) correspond to the damping by scattering processes (“Landau damping”), rjrj -H- 
XX rjx -H- TJX, respectively. The variables x and y are the cosines of the two nontrivial 
angles, i.e. x = ?/ = so that X(a, 6, c) represents integral over angles. 

The integrals cannot be performed analytically in general. For the special case p —)■ 0 
under consideration, we can, however, hnd analytic approximations, see appendix |A.2| . 


A.2 Analytic approximation for |p| = 0 and cj —> 0 

Here we derive approximate estimates of Hq (o;) for zero spatial momentum and arbitrarily 
small energy. 

In (95), if the energy po = is smaller than 2M^ + M^, then the decay term V doesn’t 
contribute as the decay process is not allowed by the energy conservation. Therefore let 
us focus on Landau damping terms (97) and (98). The integral over flk does not vanish, 
only if there exists an overlap between the integral region of ffk and the region allowing 
Oj < I (with j = 1,2) due to the step function 6^(1 — aj) from the angle integral in (176) 
(now with a = aj). 

Let us hnd the range of flk that allows a| < 1. Using the dehnition of aj given in (109) 
(with p = 0) and introducing dimensionless variables as 


XVJ.^ 


_ oj 

‘ ” A4 - ® ^ jip - ^ jiP’ “ ” 2Af|’ 


( 110 ) 


the range of t for < 1 reads 

tj_ < t < tj+ (111) 

with tj± being solutions of a| = 1, given (up to the hrst order in inhnitesimal a;) by 

ti± = {u — l)s ± \/u{i 


■u: 


t2± — 


2)(.2-i) 

^1 — M + (2m — l)s^ ± s\/u{u — 2)(s2 — 1) 

(u — l)s ± \/ u{u — 2){s‘^ — 1) 

1 + u — (2m — l)s^ ± s\/u{u — 2)(s^ — 1) 


+ 0(u^), ( 112 ) 


2m — 3 




(113) 


First of all, in order to have real tj±, it should be that u >2 (i.e. > 2M^). If m < 2 

(i.e. Mn < 2M^), then it follows that a| > 1 and integrals vanish. When m > 2, one can 
see that t 2 ± < 0 for arbitrarily small a) since (m — l)s > ■\/u{u — 2){s^ — 1), so there is no 
overlap between the integral range of flk and the range allowing < 1. S 2 thus vanishes 
due to the angle integral^ 

^^This can be physically understood in the following way: S 2 terms correspond to Landau damping 
rjx ^ VXi which effectively becomes decay x ^ VX for w —>■ 0 and this process is kinematically forbidden 
for on-shell particles. 
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For 5i corresponding to Landau damping rjr] -H- xXi which is the only potential non¬ 
vanishing contribution, one can show that M^ti± > for arbitrarily small u, so that the 
integral over ffk becomes 


dflk ^i) ^(1 — ^i) — 




dflk Fsi (flk, ^i) 




= T 


UJ ^0 rji 

= (jJT 


-log 


log 


/b(—M^ ti+)/B(—M^ti_ — f2i -|- O') 


d 


/s(—f2i_ — ffi) 


(114) 

(115) 

(116) 


./n(—f^i-)/B(—f^i+ — f^i). 

In the second and third equality we have analytically performed the integral and taken 
limit a; —)■ 0 to obtain a result up to the hrst order in u. Here 


f2i± = M^ti±{oj = 0) = {u — l)s ± u{u — 2){s‘^ — 1) 


(117) 


which depend on f2i and Mrj through s and u. 

Let us now consider the fli-integration of ( 116| ), i.e. 


T Two 

d^i ——— log 


-fli — Hi) 


f b{—^i-) f b{—^i+ — f^i) 


(118) 


This integral does not seem to allow an exact analytical result, so we try to get an approx¬ 
imate estimate of this integral. 

First, when u is close to 2 we use the fact that the integrand has a pealf^ around s = 2 
(i.e. Hi = 2M^) as shown in hgure 13 and £x the logarithm at the peak to pull it out of 
the integral. Then we can perform the integral (with lower limit Mq instead of M^) to get 


— 2a; T log 


/s(—H i+)/b(—H i_ — Hi) 

/s(—H i_)/b(—H i+ — Hi)_ 


/b(^o), 

Q.\=Mq 


(119) 


where we have put an overall factor 2 to compensate the change of the lower limit of the 
integral. Indeed numerical analysis shows that the above expression with Mq = 2M^ is a 
good approximation of the integral for 2 < M^/M^ < 3. 

On the other hand, if u becomes larger, then Hi+ and Hi_ approach the original inte¬ 
gration boundary (see hgure [l4| . Therefore for sufficiently big u, one can take Hi+ —)■ oo 
and Hi_ —)■ H||^ = max[M,j — Hi,M^] and then it is possible to analytically perform the 
integral giving 


uT 


log(l-e-^) -2 /b(M,) 


V, - 2A4 + HM, - M^] 

+2fB{M^) log 


l + /n(M, 




1 + /b(M, 


( 120 ) 


15 


To see that the integrand has a peak, note that the logarithm is zero at s = 1, at which r2i+ = tli-, 


while 


d f2i 


falls with s > 1. Furthermore, we numerically verified that the position of peak 1 < s < 2 


little changes with temperature T. 
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Figure 13: Integrand in ( 118| ) (multiplied byT) as a function of s for u = 4 (i.e. Mn ~ 
and T = lOM^. Note that the function plotted here (the vertical axis) is dimensionless. 


We numerically checked that the above result agrees with the exact integral to a factor of 
0(1) for M^/M^ > 3. 

Altogether we obtain^® an approximation of n^(a;) for inhnitesimally small u: (to a 
factor of 0(1)) 


^^o(n;) h^T 
hm —~ 


2 log 


/b(-^1+)/b(-^1--^i) 

/s( —f^l+— 




/b(2 M, 


cd^O U) 


8(27r) 


log 1 -e T - 2 /b(M^) 


Mr,-2M^ 

2T 


+ log 


fiMx) 




77 ) 


if 2M^ <M,< 3M, 


if 3M^ < M, 


and zero otherwise (i.e. if < 2M^). When My. <C T, the above expression (i.e. 

the second one) can be approximated, to leading order, as 


, mo(a;) /2T, M„\ , 

lim —~ - -log —k Q(i)_ 

uj 8(27r)3 ^ 


( 122 ) 


B The self-energy H, from setting-sun diagram 

A. 

In this appendix we calculate the spectral self-energy fl” for the self-interacting scalar par¬ 
ticle X from the setting-sun diagram shown in hgure In particular we find approximate 
estimates of fl” for the on-shell case (i.e. when the external momentum is on-shell) in 


appendix B.2 This is needed to determine T^ in ( 86 ) in section 4.2.2 For p = 0, If has 


already been calculated in [75] and here we generalize the results to arbitrary p 7 ^ 0 . 


far we have considered only 5i(a;) in (95) and have to include 5i(—w) to get 11 . Inclusion of this 
term gives rise to an overall factor 2 since 5i(w) = —5i(—w) 


(121 
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Figure 14: Integration region (shaded) in s-t plane restricted by the angle integral (i.e. af < 1) 
foru = I (i.e. ~ The dashed curve is the solution of af = 1 consisting of ti±, while 

the solid lines correspond to the original boundary of the integral region when not taking into 
account the angle integral (see figure witF flp = . 


B.l General expression 

has a general ex¬ 
pression 


The spectral self-energy If from the setting-snn diagram in fignre 


15 


nip(po) = --fj 


f d(^q d'^k dH 


-q-k- l)p^{q)p^{k)p^{l) 


X 


(27r)4 (27r)4 (27r)4 

(1 + /b(?o)) (1 +/^(fco)) (1 + /b(^o)) —/s(q'o)/b(^o)/s(^o) , (123) 


which can be derived in the same way as (94). After integrations over q and some angles 
of momenta k and 1, we get 


n„(po) = -!(»f'(po) + 5M(po) - - 5W(-po)) 


(124) 


with 


O|j'(po) 

SP(Po) 


|•po-2 M^ pU 

TT . ^ .J iPo) / -A)X(a,6,c) (125) 

24(27r)^ J 


^2 poo poo 

37r ^ .J {po) / dill / (if2k^s(f2k,fli, A)X(a,6,c). 

24(27r)° J Jq+ 


(126) 


39 
















Figure 15: Setting sun diagram for A^x^/4! interaction. Dashed line is x 


Here 


Fd(Qk, ^ 1 , A) 

Fs(f2k, A) 


(l + fsi^k)) (l + /s(fil)) (l + fsiA)) — /B(fik)/B(^l)/B(^),(127) 
(l + /s(^k)) (l + f b{^\)) f b{.A) — /B(fik)/B(^l)(l + /s(^));(128) 



= max [po — Ill ± , 

(129) 

A 

fik f^l POf 

(130) 

X(a, b, c) 

- [Xr 

7-1 7-1 ^yi{x,y,z) 

(131) 

z 

= a + by + cx, 

(132) 

l{x,y,z) 

= (1-1/^) - 

(133) 

a 

yl2_p2_]^2_j2_^2 IpI IpI 

2|k||l| ■'’=|k|’"= 111' 

(134) 


The notations here are the same as in appendix A.l The integration limits imply that 
V^\po) vanishes for po < 3M^, due to the energy-momentum conservation in the decay 
and inverse decay processes y O xXXy which are not allowed for on-shell particles. 


B.2 Approximate estimates for on-shell quasiparticles 


Here we derive approximate estimates of the self-energy ( |124 ) for the on-shell particle with 
arbitrary spatial momentum. In the on-shell case (i.e. po = hip), we hav^^ 


nip(SJp) = -3.^ 


dfli / dflkFs{flk,^h ^)F{a,b,c). 


(135) 




'Q+ 


^^Note that V terms and po) vanish for on-shell. 
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Figure 16: Division of the integral region in (137). 


Using the result on the angle integral given in (199) in appendix C.2.2, the above integral 
becomes 




dQi / dQk Fs(f2]i, Qi, A) X(a, b, c) 


'My. 


'n+ 


271 / dQi / dQkFs{Q]^,Qi, A) 
-\-2 (-—f dill |1| f dllk 


vIpI/ 


'My, 


Tl f 

IPI J My, 






(136) 


I 

II 

III (137) 


In the second line, symmetry with respect to flk -H- f2i has been taken into account, 
resulting in overall factor 2. Here I, II and III refer to the corresponding integral regions 
in plane as shown in figure 16 


The first line (integral over region I) can be analytically performed, giving 



27r / dfli / dflk-Us(f2kab ^) 


271 




2irry7rV6| if ^ 0, 

0 if ^ CO. 


(138) 

(139) 

(140) 
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When fip = 
mode, so that 


J^s and (139) corresponds to the result of the entire J^s for zero 


nxo(^x) 


—Svr 


^Xl{27^T^) 

24(27r)5 




(141) 


Now let’s turn to integrals over region II and III. Concerning II, at most one integral 
can be performed, giving 


XM = 2(^ 


PI/ J Mt, 




AttT n 

IpI J 


d^\ |1| (— f^p)) log 




The above integral and the one over region III, 


jrTI] ^ ^ 




IPI J 


dill [ " dflk jA^-M^F,{n^,nuA), 

/ A/f 1 o o. V 


' I 


(142) 

(143) 


(144) 


can not be done analytically. Therefore we try to obtain approximations of the integrals 
over II and III by considering cases of flp < T and flp > T separately, under the assumption 
< T. 


B.2.1 np<T 


- By expansion of exponential in distribution functions: 

expand the exponentials inside fs in (143)-(144), i.e. ~ 

algebra we get 


For flp < T one can 
T/Qi , and after some 


with 




G{x) 

Gi{x) 

G2{x) 




1 - t 


X^ 


dt 


dt 


L 
L 


ds 


ds 


^J{t + s — ly — x"^ 
t s {t + s — 1) 
\/{t + s — ly — x"^ 
t + s — 1 


(145) 

(146) 

(147) 

(148) 

(149) 
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Here we have introduced variables t = s = and x = M-^/Vt-p. From above 

expressions for the G functions it follows that 


G(0) = ^, G(1)=0, 

(150) 

Gi(0) = ^, Gi(l) = 0, 

(151) 

^2(0) = y g^ 2 (i) = 0. 

(152) 


We have numerically checked that for x < 1/2 (i.e. hip > 2M^) the G functions in (147) 


(149) change rather slow with x and that for x < 1/3 (i.e. flp > 3M^) they are little 


different from the values at x = 0. Thus for 3My < flp < T we can replace the G functions 


in (145)-(146) by their values at zero, giving 


O 7r2 

47rT2^-0(l), 

IpI o 


IpI 




VLr 


— Oil) + ^0(1) 
6 ^ ^ 2T ^ ' 


Together with (|139|), the entire result for the self-energy (|135|) reads 

log (l 




37r- 


^A2 (27rT2) 
24(27r)5 


Lio 


1 nil 

1 — e T 


+ 


Qr 


IPI 


(153) 

(154) 


TT 


Qr 




2T 


(155) 


Another way: The function fli, A), given in (128), has a maximum at fli = M. 


and Ilk = (also at flk = and fli = flp, by symmetry fli -H- flp). Around this point, 
Fs is dominated by the term /^(fli)^ for fli < T. Then one can use /^(fli) ~ T/Qi to see 
that |1| /^(fli)^ has a peak at |1| ~ M^. Therefore the integrand in (142) has a maximum 


My and in a similar way one can show that the integrand in (144) has a peak when 


at |1| 

replacement |1| —)• My in (142) and 
analytically performed, giving 


M2 My. Using these facts, we can obtain approximations to (|142|) and (|144|) by 


M2 ^ My in rtl44F Then the integrals can be 




dfii|i| / dn^Fs{n^,ni,A) 


I 


'On 


o f 27rMx 

V IpI 


dfli / dflk .f/A) 




[AnMyT^^ 


V 


iPl 


'vr^ (log (1 - e-^-/^)y (log (1 - e-^-/yy 

y ^ 2 2 


log 


1 _ e-^x/T 


log 


- Li2 


1 _ e-W/'T 

X — g—Op/T 


(156) 
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and 


27r 


jrpi] = — dill 


IPI JMv 


' H“ ^ p— 


dfik ^A'‘ - Ml F.(nk, Oi, A), 


(167) 


2ttM /"^p /"^p 

^ ‘ dill dnkFs{nk,ni,A), 


IPI J J 

^2 (log ,_^_Op/T J fl-e-^F/T 

“ - 2 - 


2ttM^T‘^ 

IpI 


1 — g ^x/'^\ / 1 — g(-^x Op)/T 

log ( ^- __o_ fT j log 1 ;- _o 1 - Li 2 


+Li2 


1 — (^—^pjT 

1 _ g-A^x/'T 
X — g(~A^;,.—f2p)/T 


1 — 1"^ 

X _ g~Op/r 

Lio 


X — Op)/T 


X — 6 Op/r 

g-Mx/T _ g-Op/T 
g-Mx/T _ X 

'X - g-A^x/r 

+ Li2 


X — 6 Op/T' 


(158) 


These expressions have more involved forms than (153)-(154), which can be seen as ap¬ 
proximations to the above expressions. 


B.2.2 f2p>T 


In this case the integrands in (142) and (144) are maximal at |1| rsj T and J^ T, 
respectively. Therefore we can obtain approximations by replacing |1| —T in (142) and 


A? — —)■ T in (144). Then we can analytically perform the integrals to obtain 


27r 


j-m = 2 — / dfli|l| / 


IPI 

27rT 


'Mv 


'O 


■(iPi 


P 

oo 


/ dill / dflkT;(f^k,f^i,4l) 

'Mq Jup 

TT^ (log (l — (log (l — 


log 


1 — e ^p/"^ 
X — q-MqIt 


log (e 


,-MolT _ -Oj 


/^) - Li2 


X _ g-A^o/T 

X — g—Op/T 


(159) 
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and 


27r 


= A / dn, 


IPI 7m, 
2ttT 


' -|-0p—j 

/■f^n 




(160) 


7f2i / dUk T^s (^k) ^1) ^) 

Pl J Mq J Mo-\-^p — ^i 

7r2 (log i_e-^p/^ J /I - 

y+ ^- 4 -^-(^P + ^o)log' 


27rT3 

IpI 

+ log 
+Li2 


1 _ e-Mo/T\ ^ _ g(Mo-np)/r 

j log 


1 — g—f^p/T 

1 _ g-A^o/T 
X — g(—A^o—Op)/T 


— Li*; 


X — g( Afo Op)/T 

1 — e“Op/^ 

' X - g(-Afo-np)/T 


1 — 6 Op/T' 

,-Mo/T _ g-Op/T 
g-Mo/T _ X 

X _ g-ATo/T 


+ Li 


1 — 6 Op/T' 


(161) 

Here we have replaced in the low limits of the integral by some constant Mq, which we 
choose Mq ~ T/2. For flp S> T, from above we get 


- 2 vrT^ ( 4 - log(e 


^0 , / Mo/T _ 


1 ) 


(162) 

~ 2tiT‘^0{1), (163) 

and this is the result of the whole Fg since ~ 0 for flp S> T (see (140)). 

To sum up, when -C T the spectral self-energy for on-shell y-particle with any 
momentum p can be approximated as 

i\l {271T^) 


n-p(^^p) 


-Svr- 


0 ( 1 ). 


(164) 


24(27r)5 

C The angle integral in setting-sun diagrams 

In this appendix we estimate the angle integrals for setting-sun diagrams I{a,b,c), see 


(106) and (131). The angle integral has a general form 


X{a,h,c) = 

/‘dx 

7-1 7-1 ^yi{x,y,z) 

(165) 

I{x,y,z) = 

(l-x^){l-y^)-{z-xyy, 

(166) 

z = 

a + by + cx, 

(167) 

a = 

H 2 -p 2 -k 2 -i 2 -M 2 ^ IpI IpI 

(168) 

2 |k||i| |k|’^ |ir 


Here H is a function of Pq, flk and fli, see e.g. (130). First we consider the cases with zero 


external momentum p = 0 (on-shell or off-shell) and then turn to the on-shell cases with 
non-zero momentum p 7 ^ 0 . 
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C.l Zero mode (p = 0) 

In this case we have b = 0 = c and z = a, so that 

(169) 

If < 1 and < 1, there exists —1 < ?/± < 1 that solve I{x,y) = 0, so that one can 
write (see hgnre 


= -{y - y-){y - y+) 


with 


y± = ax ± a/(1 — x‘^){l — a^) 


(170) 


(171) 


One can easily see that |j/±| < 1 for < 1 and < 1. Then assnming < 1 the angle 
integral can be performed, giving 


X 


-1 


dx 


e{I{x,y)) 


L 

/: 

2ti. 


dx 


-1 ^ \/i{x,y) 

ry+ X 


dy 


'y- 


\/-{y - y-){y-y+) 


dx 71 


(172) 

(173) 

(174) 

(175) 


Note that in the second line the ^/-integral gives tt irrespectivly of x, thongh y± depends 
on X. If > 1, there exists no real y± for < 1. This implies I{x, y) < 0 for < 1 and 
y"^ < 1, giving vanishing resnit for the integral due to d{I{x,y)) in the integrand. 

The hnal result for zero mode is thus 


I = 2776(1-a^). 


(176) 


Note that a depends on po, |k|, |1| and by (168), so that the step-function 9(1 — a?) in 


(176) restricts the integral region on Ok-Oi plane to the one allowing < 1. 


C.2 Non-zero mode (p ^ 0) 

Here we generalize the zero-mode result obtained above to the cases of non-zero mode. 
First we consider generic cases, which include off-shell as well as on-shell, and then focus 
on the latter, which is needed to evaluate T^ in (86) in section (4.2.2). 
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C.2.1 Generic cases 

For Ip I 7 ^ 0, 2 ; is given hj z = a + by + cx with 6 7 ^ 0 7 ^ c and /(x, y, z) can be written as 
Hx,y,z) = -{l + b'^ -2bx){y-y+){y-y-), (177) 


where 


y± 


{x — b){a + cx) ± -^(1 — x 2 )(l + 6 ^ — 2bx — (a + cx^) 

1 + b"^ — 2b X 


(178) 


Let us find the conditions, under which there exists — 1 < |/± < 1 assuming |x| < 1. The 
conditions that y± be real and that it satisfy — 1 < |/± < 1 require 

{a + cxY < (1 + 6^ —26a;), (179) 

{x — b)'^ {a + c xY < {l + b‘^ — 2bxY, (180) 

respectively. The second condition above is obtainecj^by requiring — 1 < < 1. Since 

0 < (x — 6)^ < 6^ — 26x + 1 for |x| < 1, the second condition also follows from the first one 

with |x| < 1. Thus the condition for — 1 < < 1 is 


(a + cx)^ < (1 + 6^ — 26x) and |x| < 1. (181) 

The range of x satisfying this condition is 

X = (-1,1) n (x“,x+) (182) 

^®Since I{y = ±1) < 0, j/_ and y+ must lie in < — 1 } or {—1 < y-,y+ < 1} or {1 < y_,?/_|_}, see 

figure [T7| Among these possibilities, only when —1 < < 1, the integral does not vanish and only 

this case fulfills — 1 < (?/_ + y+)/2 < 1 . 
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where are solutions of {a + cxY = {1 + — 2b x), 


x^ = 


and they are real when 


-b — ac± 1/52 + 2abc + + b‘^c^ 


6^ + 2abc + + b‘^c^ > 0. 


Furthermore it follows from (168) that 




b + 2abc + c + be = 


iPl 


(183) 

(184) 

(185) 


The right hand side is positive over the integral region in (125) and (126) since it has been 
chosen in such a way that holds, by using for the integral limits. Thus x^ 

are always real on the integral region. 

If X = 0, then the angle integral vanishes. If X 7 ^ 0, we define x_ and as a lower 
and upper bound of X, respectively and then the angle integral can be performed as 


X{a,b,c) = 


1-1 7-1 \/I{x,y) 

f dx /■^+ 


dy 


lx Al + b^-2bx Jy_ ^y-{y-y+){y-y-) 
dx 


r-a;+ 


Jx- a/1 + 5^ — 26x 

TT ^a/1 + 6^ — 2bx- — a/1 + 5^ — 2bx+'^ 


(186) 

(187) 

(188) 

(189) 


Here X- = max[—l,x ] and = min[l,x’''] provided that X 7 ^ 0. Since x^ are functions 
of a, b and c, the condition X 7 ^ 0 can be reduced in terms of these parameters, which are 


functions of energy, momenta and masses of particles through (168). In the following we 


consider on-shell cases and obtain analytical results for the angle integral. 

C.2.2 On-shell cases 

As a specific example, we apply the results obtained in the previous subsection to the 


diagram given in hgure 15, when the external four-momentum is on-shell (i.e. po = = 


-|- M^). By introducing 7 = 1 ^ the parameter a can be expressed as (using (168) 
and (130)) 


a = 'j'^bc + a /(1 -I- 7262)(1 -|- y^c^) — ^1 -|- 7 ^ + 7^^^ + &\/l + , (190) 

which has properties 


a{b = 1 ) = —c, a{c = 1 ) = —b. 


(191) 
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(192) 


Using the second property above in (183), we get 


and one can show that 


Therefore 


x^(c = 1) = ±1 

|a;^| >1 if c < 1, 
Ix^l <1 if c > 1. 


X = (—1,1) n (a: , x"*") = (x_, x+) = 



if c < 1, 
if c > 1. 


Then in (189) we have 




a + cx 


±1 


if c < 1, 
if c > 1. 


(193) 

(194) 


(195) 


(196) 


In the second line on the right hand side we have used the fact that solve l + b'^ — 2b x = 
{a + cxY- Noting that a + cx~ < 0 and using properties for c > 1, 


a + cx 


<0 if 6 > 1 

= 0 if 6 = 1 

>0 if 6 < 1, 


(197) 


we obtain the hnal expression for the angle integral (189), 


TT 


X = 


^a/ 1 + 5^ — 26x — 


= 


27r 

27Tb-^ 
271 c~^ 


if 6 < 1 and c < 1 
if 6 > 1 and c < 1 
if 6 < 1 and c > 1 


(198) 


27r I — if 5 > 1 and c > 1. 
IpI 


In the last line we have used (185) and assumecf^ A > My. (i.e. flk + + My. in 

(126)). In terms of energies and momenta, the result above reads 


X = 


27r 


if 

if flk f^p ^ 
if < flp < flk 

if flk, < f^p and < Ilk + 
otherwise. 


IpI 

IpI 


(199) 


Qr 


0 


As mentioned before, A > holds over the integral region on flk-Hi plane in (126) by imposing the 
integral limit 
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